2, 2008 8:15 WSPC/INSTRUCTION FILE cflslet 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



Color-flavor locked strangelets in a quark mass density- dependent 

model 



X. J. Wen, 1 - 2 ' 3 * G. X. Peng, 2 > x P. N. Shen 1 ' 2 

1 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China 
2 China Center of Advanced Science and Technology (World Lab.), Beijing 100080, China 
3 Graduate University of Chinese Academy of Sciences, Beijing 100049, China 



The color-flavor locked (CFL) phase of strangelets is investigated in a quark mass density- 
dependent model. Parameters are determined by stability arguments. It is concluded that 
three solutions to the system equations can be found, corresponding, respectively, to 
positively charged, negatively charged, and nearly neutral CFL strangelets. The charge 
to baryon number of the positively charged strangelets is smaller than the previous result, 
while the charge of the negatively charged strangelets is nearly proportional in magnitude 
to the cubic-root of the baryon number. However, the positively charged strangelets are 
more stable compared to the other two solutions. 
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1. Introduction 

At high densities, quark matter, or its finite lumps, the so-called strang elets , may 
form the color superconductor by chiral-symmetry violating condensate!!^ Most 
probably strange quark matter (SQM) is more stable than the hadronic mat tersP 
which may have far-reaching consequences astronomically and c osmologically.™ ^ 
It has potential applications to the astrophysics of compact starsj ^^ l lO l ll l l^ l anc j 
to high energy heavy ion collisions ff^ Recently, the prop erties of CFL phase have 
been investigated by many authors with various models j l4 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 
Rajagopal and Wilczek demonstrated that the CFL quark matter is automatically 
neutralized without any requirements on electrons!^ Madsen studied the positively 
charged CFL strangelets, and found that the CFL stran gelets m ay be more stable 
than the ordinary strangelets without color-flavor locki ng R^^ Alford investigated 
the new gapless CFL phase in neutral cold quark matter P'^'^ Lugones and Horvath 
derived an approximate equation of state (EOS) to lowest order in A and m and 
found the effects of pairing on EOS and the window of stability for CFL strange 
matter P3 
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In the standard MIT bag model, quark masses are constant. As is well known, 
however, the quark mass changes with environment, or, in other words, it depends 
on density. Such masses are usually called effective masses. | 28 | 29 | 30 | 31 | 32 | Effective 
masses and effective bag constants for quark ma tter have been extensively dis- 
cussed within the Nambu-Jona-Lasinio modelP^ In recent years, the quark mass 
density-depe ndent model (QMDD) has been su ccessful in describing the unpaired 
quark matter! a4 l a5 l ^ l a7 l ^ l ^ l 4 » l 41 l 42 l ^ l » Very recently, we also studied CFL 
strangelets within the framework of bag modelPSl Now in this paper, we study how 
the density dependence of quark masses influences the properties of CFL strangelets. 
It is found that for a given baryon number, there are also three kinds of CFL 
strangelets which are, respectively, charge-positive (slet-1), negative (slet-2), and 
nearly neutral (slet-3) . These multiple solutions are determined by the multiple val- 
ues of quarks Fermi momenta satisfying the self-consistent condition. The finite-size 
effects of strangelets with fixed Fermi momenta influence the fractions of quarks as 
well as the charge of strangelets. For slet-1, the ratio of the charge to squared cubic- 
root baryon number is only half of that in the pure bag model, while the charge of 
slet-2 is proportional, in magnitude, to the cubic-root of baryon number. For the 
same parameters, the slet-1 is more stable. 

This paper is organized as follow. In the subsequent Sec. [2j we briefly give the 
thermodynamic formulas used for the study of CFL strangelets in the mass density- 
dependent model. Then we present the numerical results and discussions in Sec. [3] 
The Sec. H] is a short summary. 

2. Thermodynamic formulas and mass density dependence of 
quark masses 

In this paper, the CFL phase consists of up, down, and strange quarks in three-color 
QCD. The symmetrical color-flavor locked phase needs an attractive interaction 
between two quarks near the Fermi surface with equal and opposite momentum. The 
Fermi surface is vividly fixed in the momentum space. The condensate of Cooper 
pairs, associated with the group of ud, us and ds, spontaneously breaks the color 
gauge symmetry. At high densities, quarks of three colors and three flavors are 
allowed to pair and hav e th e same Fermi momentum. 
In the pairing ansat J^31 

(CC75^> ~ A 1 e a ^ 1 e abl + A 2 e a f 32 e ab2 + A 3 e a ^e ab3 , (1) 

we consider the case A3 A2 = Ai = A for common CFL phase. Here ipZ IS 
a quark field with color a — (r,g,b) and flavor a — (u,d,s). T he therm odynamic 
potential of the color-flavor locked strangelets can be written as H5123147] 

O = ^ Q[ t i + fipair, v + B , (2) 

i 

where the bag constant B is appended as normally done. Eq. (|2|) is derived from 
microscopic models under the condition that A//i is small!^ For pairing contri- 
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bution, we include the volume term fi pa ir,v ~ — 3A 2 /i 2 /7r 2 , which will appear in 
the number density, energy and pressure etc. The finite-size contribution to fi pa ; r 
is assumed small and reasonably neglected!^ The quantity ft = + fid + /i s )/3. 
Q{ t i denotes the normal quark contribution from flavor type i ( i = u, d, s ), i.e., 



fi f)i = -T hxi 1 ± exp 



T 



n'tip, rrii 1 R)dp . 



(3) 



The integral upper limit pf me ans Fermi momentum, which is equal for u, d, s 
quarks in CFL quark matter!^ When T = 0, this expression can be simplified, 
giving 

i-PF 

fif,i = / (si - fJ.i)n' l (p, m t ,R)dp, (4) 



where the density of states for strangelets from the multi-reflection theory-^ is 

P 



(5) 



Here V — jirR, S = 47ri? , and C — 8irR. The function /? and / 2 G are respectively 
referred as surface term I3ESI and curvature term I 50 | 51 | an( j gj ven by 





1 


- 

rriij 






1 


rriij 


12tt 2 



-1 ( p_ 

mi 



1_ JP-ct-if A 
2m, 



(0) 
(7) 



For convenience, we define a parameter 0i = arctan(pp/™i). If the ther- 
modynamic potential density in Eq. (JJ) is divided into three parts i.e., fif j = 
fi^ + 4 fif + "/fr^f , then the volume term Qj , the surface term fif, and the curve 
term fif can have analytical forms as, 



I67H ( 3 mi 



- In [tan^ + sec fa] 
— tan fa sec 3 — tan 3 fa sec <fo | , 

3 ( 

fif = -|t^| \ In Nn + sec fa 



3pi T.7T 2 

(—-(%) tan 0i+tan0i- 

m, L 2 



+ : 



■ sec fa [tan fa + sec </>i(7r — 2 fa)] — 
fif = | In [tan & + sec fa 



3/Xj 



7T 

- tan sec fa — 2 sec 3 ^ ( — — ^ ) + 



1 7T 

- tan ^ - tan 2 fa(- - fa) + < 

17 



7T 



(8) 



(9) 



(10) 
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The common fermi momentum p-p is a fictional intermediate parameter, which 
can be determined by minimizing the thermodynamic potential f2, i.e., 

on 



or, explicitly, 



dpi 



= o. 



(ii) 



p\ + mf - [ii 



= 0. 



(12) 



By differentiating fltA with respect to the chemical potential //$ we obtain 



nf,i = 3ip|/(67r 2 



8tt 2 R 



2 r 



tan < 



] tan z 



87T 2 i? 2 

Hence the number density for flavor i is 

rii = UfA 



1 7T 

-tan^ - (- - (^)tan 2 . 



(13) 



(14) 



The second term is from the paring effect. We should accordingly note that the 
density depends directly on the chemical potential jj, and the paring parameter A, 
not merely the Fermi momentum "pf" ■ The corresponding energy density can then 
be expressed as 



E = ft + [liUi 

i 

= X [fif,j + Hi(nfA + n pair! y)] + ilpair^y + B . 



(15) 



In the conventional bag model, quark masses are constant. As is well known, how- 
ever, quark masses vary with environment. In fact, not o nly quark masses will change 
but also the coupling constant will run in the medium.^. Effective masses and ef- 
fective bag constants for quark mat ter have been extensively discussed, e.g., within 

looi iroi 

the Nambu-Jona-Lasinio model 5 ^ and within a quasiparticle model. — Therefore, 
in recent years, the quark mass den sity-dependent model has been show n to be suc- 
cessful in the study of quark matter. IM13 5 | 36 | 37 | 38 | 39 | 4() | 41 | 42 | 43 | 44 ]45l The question 

now is how to parameterize the density dependence of quark masses. In principle, it 
should be treated dynamically and be consistent with th e overall chir al symmetries 
of QCD. In studying the unpaired phase of quark matter, [ 38 | 39 | 40 | -| 5] equivalent 
quark mass is 

D . . 

m q0 + . (16) 



m„ 



The baryon number A is connected to the quark number N u , Nd and N s by 



A=^(N U +N d + N s ) : 



(17) 



2, 2008 8:15 WSPC/INSTRUCTION FILE cflslet 



Color-flavor locked strangelets in a quark mass density- dependent model 5 

and the chemical equilibrium required 

Md = Ms = Mu + Me ■ (18) 

When the baryon number j4 -C 10 7 , electrons can't exist in a strangelet because 
the electron Compton wavelength exceeds the sphere radius R. So the electron 
"effective" chemical potential is zero in strangelets, which is consistent with the 
viewpoints of Ref . IS"4l and 1551 to some extent. Figure 11.1 of Ref. 54 and Figure 3 of 
Ref. [55]point that when the mass of quark matter is less than 10 9 GeV, the electron 
could be outside the quark bag. So when /i e — 0, Eq. (|18p becomes fi u = jid = Ms- 
With this relation, the chemical potential can be obtained from Eq. |[T2|) for a fixed 
A and a given radius R. 

The radius of a strangelet with a given baryon number is determined by mini- 
mizing the thermodynamic potential with respect to the radius, 

M = °- ^ 

Or, equivalently, by setting the pressure 

P- Q -— dmi 80 B (20) 
3 dR Uh drib drrii 

to zero, i.e., P = 0. The first two terms on the right hand side of Eq. (|20|) is 
normal, while the third term is special when quark masses are density-dependent. 
As has been shown in literature, this extra term is necessary to satisfy the basic 
thermodynamic require ment, i.e., the energy minimum must be exactly located at 
the zero pressureP^H 

According to Eqs. (JH]), (J9j> , and (fT0|) . the derivative of thermodynamic potential 
density with respect to mass can be divided into three parts as, 



(21) 



9to~ = An 2 X tan sec ^ _ hi(tan</>i + sec fa)j, 

= — ^i m i j tan ^ _ gec ^ _ _ sec 0. _|_ i n (tan fa + sec fa ) 
ami 87H ^ 

- (1-860^)^ + 2— [fa- tan fa] \ , (22) 
rrn J 

dn c gi-nii ( /3//i 

— — = -— r < 4 ln(tan fa + sec fa) + 2 tan fa 

ami 9D7r z ^ V m t 

+ 7r I 2 tan 2 fa sec fa — 4 sec 

- 2 fa ( 2 tan 2 fa sec fa + ^ - 4 sec 0, - 3 ^ tan2 M 1 . ( 23 ) 

Because no electrons are included in the system, finite-size CFL strangelets are 
charged. The net electric charge is Z = |n u — ^rid — ^n s in unit of the electric 



— 2 sec < 
tan 2 fa 
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charge of an electron. In principle, the Coulomb effect should be included, though 
it is small compared with the strong interaction. We include it in the numerical 
calculations. 

For the charged CFL strangelets, the Coulomb energy is, 

1 aZ v laZ 2 

Ec ° ul - tq— + > (24) 

where the fine structure constant is a ~ 1/137 and Zy is the volume term of the 
total electric charge Z of the CFL strangelets. 



3. Numerical results and discussions 

We now discuss the parameters adopted in our calculations. The current quark 
masses are taken to be m u0 = 5 MeV, m^o — 10 MeV m s0 = 120 MeV except for 
specific indications somewhere. The super-cond ucting ga p varies from several tens 
to several hundreds of MeV in previous papersP^E^B] j n ^his p a p 6 r we take the 
value A = 100 MeV, as in Ref. H3 
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Fig. 1. The stability window of D 1 / 2 and 
B 1 / 4 . We take m u0 = 5 MeV and m d0 = 10 
MeV. The lower bound is given by the full line, 
while the upper bound is given by the dashed 
line for m B o = 80 MeV and by the dotted line 
for m s0 = 120 MeV. 
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Fig. 2. Energy per baryon of the positively 
charged CFL strangelets for A = 200. The 
solid, dashed, and dotted curves stand for 
D 1 / 2 = 80 MeV, 40 MeV and 0, respectively. 
The bag constant is B 1 / 4 = 150 MeV. Fermi 
momenta are also indicated on the right axis. 



It is known that the lowest energy per baryon in ordinary nuclei is 930 MeV 
for iron. In order not to contradict with standard nuclear physics, the energy per 
baryon at zero temperature should be greater than 930 MeV for two flavor quark 
matter, and less than 930 MeV for three flavor quark matter so that SQM can 
have a chance to be absolutely stable. We show the stability window in Fig. [TJ The 
horizontal axis, namely D = and constant quark mass m q = m q o (q — u,d,s), 
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Fig. 3. Energy per baryon and electric charge 
per A 2 / 3 as functions of baryon number A 
for CFL strangelets (slet-1) are shown. Here 
Fermi momentum pp is in the usual range of 
250 ~ 300 MeV. The ordinary strangelets are 
also plotted. 



Fig. 4. Charge of strangelets for slet- 
2 (dashed lines) and slet-3 (dotted lines). 
The slet-3 is nearly neutral, while negatively 
charged slet-2 is proportional to A 1 / 3 , the cor- 
responding Fermi momenta are also indicated 
on the right axis. 



shows a range of (144.297, 157.3634) MeV for the constant B 1 / 4 in bag model. If we 
require the energy per baryon less than the mass of nucleons E/rih = 939 MeV, we 
can derive the same (meta)stable range of B as the result in Ref.|31 The vertical axis 
shows the previous range of (154.8278, 156.1655) MeV for \[D in QMDD model.ES 
The parameter pairs {D 1 / 2 , B 1 / 4 ) under the solid line or above the dashed line will 
not be adopted for stability purpose. If B 4 ! 4 = 150 MeV, a range for D can be 
found. For baryon number A=200, the three lines of energy per baryon are shown 
respectively for _D 1 / 2 =80, 40 and in Fig. [2] The full dots denote the zero pressure 
points, which are exactly the minima of the energy per baryon. The corresponding 
Fermi momenta, indicated on the right axis, are also affected by the mass parameter 
D. 

We choose the parameters B 1 / 4 = 144 MeV, D 1 / 2 = 120 MeV in calculating the 
data for Figs. El |4] and [5] The energy per baryon and electric charge per A 2 / 3 as 
a function of baryon number A for CFL strangelets is given in Fig. [3] The strange 
current mass m s o is taken to be 80, 120 MeV respectively. We find the net electric 
charge Z w 0.15A 2 / 3 , which is half of the result in the pure bag modelP^ l n the 
same figure, the energy of the neutral ordinary strangelets is also plotted by a dash- 
dotted line. When we solve Eq. (fl2|) we find that there are three different solutions 
for the common Fermi momentum. In Fig. |4] the charge of the other two solutions 
are given. The charge of CFL slet-2 (dashed lines) is proportional to cubic-root of 
baryon number, while CFL slet-3 (dotted lines) is nearly neutral. The corresponding 
Fermi momenta are also indicated on the right axis. Specially, given A — 200 we can 
get the smaller radius R = 5.55 with p-p w 267 MeV, which is marked by slet-1 in 
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Fig. 5. The energy per baryon of CFL strangelets as functions of the radius. The mechanically 
stable radius, marked by full dots and given numerically by i? s tablei is located at the minimum, 
where the pressure is exactly zero. For a fixed baryon number (A = 200 for the figure), there are 
three solutions: Figure (a), (b), (c) respectively for slet-1, slet-2, and slet-3. 



Fig-OJa). The net positive charge Z satisfies Z/A 2 / 3 w 0.15 as mentioned above. The 
corresponding high p-p is the usual case previously obtained by other authors. With 
the Fermi momentum decreasing we can find the other negatively charged solutions 
marked by slet-2 and almost neutralized solution marked by slet-3, respectively, 
in Fig. [5]Jb) and (c). They have larger radius 7.8 fm. The corresponding Fermi 
momenta (42 MeV and 12 MeV or so) are smaller than that for slet-1 which will be 
interpreted in the next paragraph. It can been claimed that slet-1 is more stable at 
the same parameters. In Fig.[5][c) for slet-3, the charge to baryon number ration is 
Zj A w —0.0006, close to neutral state because n u w rid ~ n s is nearly satisfied. 

There exist in literature two treatments concerning the deconfined quark phase. 
One is the bulk quark matter which is infinite in volume and without regard to 
boundary of the sphere. The other is finite-size lump quark matter, i.e. the so-called 
strangelets. The ordinary quark matter can support electrons within it and is charge 
neutral in /3 equilibrium. CFL quark matter is automatically electrically neutral 
without participation of electrons because BCS-pairing requires the Fermi momenta 
of different flavor quarks to be equal. For ordinary strangelets or CFL strangelets, 
we should fall back on how the density of states depend on the quark numbers. The 
surface term has a negative contribution to the total number of massive quarks no 
matter how large the paring Fermi momentum is. The curvature term, however, 
has a positive or negative contribution, depending on the values of quark mass 
and Fermi momentum. So ordinary strangelets is positively charged because the 
number of heavier, negatively charged strange quarks is suppressed compared to 
the number of lighter up and down quarks. Similarly, the CFL strangelets (slet-1) 
is positively charged as in Madsen's work. The Fermi momenta of three cases 
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satisfy (slet-3)<(slet-2)<(slet-l). In the case of smaller pp, the surface and curvature 
terms make the number of strange quarks greater than, or nearly equal to, that of 
u/d quarks. As a consequence, the slet-2 is negatively charged. For slet-3, the very 
small Fermi momentum make the surface and curvature effect so weak that slet-3 
is nearly neutral. Although there are multiple solutions for a fixed baryon number, 
it is necessary to declare that the positively charged slet-1 is more stable than and 
the other two if the density dependence of the pairing parameters is not considered. 

It should be emphasized that the common Fermi momentum "pf" is only a 
fictional intermediate parameter in CFL matter. Different from the usual phase, the 
density of CFL strangelets depends not only indirectly on the chemical potential 
through pp, but also directly on the chemical potential. As a consequence the density 
can be large even if pp is small as long as the chemical potential is big. Therefore, the 
density is still higher than, or at most, near the normal nuclear saturation density, 
because the chemical potential is big, though pp is small. In fact, the corresponding 
chemical potentials p, of slet-2 and slet-3 are larger than that of slet-1. The results 
satisfy the validity condition that the ratio A/// is small. So the densities of new CFL 
strangelets "slet-2" and "slet-3" are still in reasonable range. We can still distinguish 
them through the difference of the electric charge Z per A 1 / 3 and chemical potentials 
fl in Fig. [6l The smaller the baryon number, the more visible the difference is. The 
Coulomb energy of the positively charged slet-1 has been shown in Fig. [7] With 
increasing baryon number A, the trend becomes placid. The effect on the free energy 
per baryon will not exceed 0.02 MeV. For slet-2 and slet-3 it will be much smaller 
than 0.001 MeV. 

It should be noted that the new solutions, slet-2 and 3, have an unusual common 



0.020 




Fig. 6. The charge to cubic baryon number 
(solid lines) and chemical potentials (dashed 
lines) vs. the baryon number for slet-2 and slct- 
3. 



Fig. 7. The Coulcomb energy per baryon vs. 
baryon number. Its effect on the free energy 
per baryon even will not exceed 0.02 MeV in 
our baryon number range. 
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momentum. When two paired quarks have a very small momentum, their global be- 
havior may looks like a boson. Therefore, the new solution may indicate that boson 
condensation or diquark condensate appears to some extent, and the formation 
mechanism needs further investigations in the future. 

4. Summary 

We have studied CFL strangelets within the framework of the quark mass density- 
dependent model. We have add the Coulomb interaction to the charged strangelets. 
It is found that the positive net charge is Z/A 2 ^ 3 w 0.15, nearly half of the previous 
result in the pure bag model. Importantly, with decreasing Fermi momentum, we 
can find the other two solutions which have larger radius than the ordinary solution 
of CFL phase. The new solutions are slightly negatively charged or nearly neutral. 
Although with small Fermi momentum, the strangelets have large chemical poten- 
tial. Due to the pairing effect, they have a comparable density with dense matter. 
The charge to baryon number of the positively charged strangelets is smaller than 
previously found, while that of the negatively charged strangelets is nearly pro- 
portional in magnitude to the cubic-root of the baryon number. However, present 
results depend on the parameter choice, and so, further studies are needed. 
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